Critically damped quantum search 
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Although measurement and unitary processes can accomplish any quantum evolution in principle, 
thinking in terms of dissipation and damping can be powerful. We propose a modification of Grover's 
algorithm in which the idea of damping plays a natural role. Remarkably, we have found that there 
is a critical damping value that divides between the quantum 0{-\/N) and classical 0{N) search 
regimes. In addition, by allowing the damping to vary in a fashion we describe, one obtains a fixed- 
point quantum search algorithm in which ignorance of the number of targets increases the number 
of oracle queries only by a factor of 1.5. 



Dissipation has generally been regarded as a destruc- 
tive foe in the arena of quantum mechanics, ruining quan- 
tum effects and greatly complicating theory. This is par- 
ticularly true in the case of quantum information science, 
which depends on maintaining delicate quantum super- 
positions and entanglement. On the other hand, dissipa- 
tion has many productive uses in the classical regime, 
from automobile shocks to toaster ovens. One there- 
fore wonders whether dissipation might be employed con- 
structively in the quantum information context [![. 

We propose here a natural application of dissipation 
in the quantum search algorithm. This algorithm, due to 
Grover j2|, is a mainstay of quantum information science. 
Given an unsorted database of N items of which M are 
target items, the quantum search locates one of the target 
items with 0{\J N /M) queries of an oracle. A classical 
search would require 0{N/M) queries. A well known 
vulnerability of the quantum search is the need for prior 
knowledge of the value of M [3, 0] ■ In the absence of such 
knowledge, the quantum search is not robust, produc- 
ing results that oscillate between target and non-target 
items. In this Letter, we show that introducing dissipa- 
tion into the search algorithm can damp out these oscil- 
lations. Strikingly, we find that a critical damping value 
emerges from the theory that divides between the quan- 
tum regime (low dissipation, , 0{^J N /M) queries) and 
the classical regime (high dissipation, , 0{N/M) queries). 

Although recently other fixed-point quantum searches 
have been developed 0, @, 0] , they are not designed to 
preserve the signature 0{yj N /M) behavior of the quan- 
tum search. When the damping is chosen appropriately, 
the dissipative approach maintains the 0{^J N/M) quan- 
tum behavior. Furthermore, when the damping is al- 
lowed to vary in a manner we describe, ignorance of M 
costs only a factor of 1.5 in the number of oracle calls, 
which is a very low overhead compared to other ways 
of handling ignorance of M. Overall, our results con- 
vincingly demonstrate the productive use of dissipation 
in quantum algorithms and provide an example where 
the appropriate amount of dissipation emerges explicitly 
from the theory. 

We briefly review Grover's quantum search Given 
an unsorted database of N items, one is charged with 



the task of finding any one of M target items dispersed 
throughout the database. The tool to probe the database 
is an oracle that indicates whether a given item is a tar- 
get or not. The algorithm begins by placing the state 
of the system into an equal superposition of all database 



states Es^i,...,^ |s > /VN = cos0/2\a) + sm9/2\P). 
Here, we have defined \a) to be an equal superposition of 
all M — N non-target states, to be an equal superpo- 
sition of all M target states, and sin 0/2 = a/ M/N. It is 
convenient to define Pauli operators in the 2-dimensional 
Hilbert space of |a) and |/3) by AT = (a| + |a) (^|, 
Y = i 1/3) (a| - i \a) and Z = \a) {a\ ~ |/3) {f3\. In 
terms of these operators, the oracle is just Z . Grover de- 
veloped a series of gates, consisting of one call to the 
oracle Z followed by one "inversion about the mean" 
E, that produces the rotation G = EZ = exp{—i9Y). 
After R applications of G, the system has rotated to 
cos(2i? -I- l)6'/2 \a) + sin(2i? -I- l)6'/2 |/3). If d is not 
too big (M <C N), and one chooses R to be near 
R{M) = [arccos(v^M/iV)/6'], then the final state will 
be close to the target \(3). Clearly, if one is ignorant of 
M and keeps applying G past R ~ R{M), the system 
will rotate past |/3). 

In terms of the system density matrix, the evolution 
p' = GpG^ implies 
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satisfies 



where the initial density matrix p — 
[Tr(pA) Tv{pZ) Tv{p)f = [sin 61 cos 6* 1]"^. 

So far, we have simply reviewed Grover's algorithm, 
pointing out how it assumes prior knowledge of the num- 
ber of targets M . Suppose that we do not know M . How 
can we find a target reliably? We propose introducing 
damping to slow the rotation in the |a), \(3) plane as the 
system approaches the target state. 

We append an external spin to the system, taking 
IV') to IV') li)- Let the external spin Pauh operators be 
Sx,Sy,Sz. This external spin will serve to indicate the 
proximity of the system to the target state. We replace 



2 



the Grover rotation G with 



U 



G 



1 + 5, 



-i4>Sy 



l-Z 1+Z 



(2) 

First, [e-"l'^y{l - Z)/2 + (1 + Z)/2] calls the oracle and 
flips the external spin if |^) has reached the target state 
{Z = -1). Second, [G(l - S^)/2 + (1 + S^)/2] applies a 
Grover rotation only if the external spin has not flipped 
{Sz = —1). Thus, the external spin stifles the Grover 
rotation as the target is approached, effectively damping 
the system. Although it appears at first that U requires 



2 controlled oracle calls, in fact it can be written using 
one controlled oracle call as 



U 



E 



i4)Sy /2 



i4>Sy/2 



Assume we follow application of U with measurement 
of the external spin, and repeat these two steps until 
the external spin has flipped. The density matrix of the 
system will have the form p ||) (|| until the external spin 
flips. The iteration produces the following effect on p: 
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Note that, since there is some probability at each itera- 
tion that the external spin will flip, Tr(p) can decrease. 

The value of <f> determines the amount of damping of 
the iteration. When = 0, there is no damping, and 
we recover Grover's quantum search. When (p — 7r/2, 
the damping is strongest - the external spin essentially 
acts as a pointer for a full measurement of Z, the system 
undergoes collapse at each oracle call, and the search is 
nearly classical. This is the limit discussed in 0, 0] ■ To 
choose the optimal value of for our purposes, we con- 
sider the eigenvalues of the matrix in (l3|). As Fig. 1 
shows, for small 0, when U is nearly G = exp(— i6'y) and 
([3]) is nearly ([1]), there are two complex conjugate eigen- 
values ~ exp(±i26') associated with the rotation and one 
eigenvalue ~ 1 that keeps Tr(/9) nearly constant. When 
(j) is large, there are three real eigenvalues. In the extreme 
case 4) — tt/2 (not shown in Fig. 1), two eigenvalues are 
and one equals cos^ 9. The eigenvalue associated with 
Tr(pX) is since the measurement destroys the coher- 
ence. The eigenvalue associated with Tr(/9(1 — Z)/2) is 
since p is associated with the non-target part of the 
system for which Z = 1. The remaining eigenvalue 
is cos^ 9 because a projection to the non-target state 
Tr(/9(1 -I- Z)/2) followed by a Grover rotation through 
9 returns a non-target state with probability cos^ 9. 

Remarkably, there is a critical damping defined by 
COS0 = (1 — sin6')/(l -I- sin0) at which all three eigen- 
values are cos <j). At this critical damping, all three com- 
ponents Tr(pAr), Tr(pZ), and Tr(p) appearing in the 
map ([3]) tend to be suppressed since all the eigenvalues 
have magnitude under 1. Thus, 1 — Tr(p), the probabil- 
ity that the target state has been found and the external 
spin has flipped, tends to increase. Fig. 2 exhibits the 
special character of the critical damping value. It de- 
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FIG. 1: Real and imaginary parts of the 3 eigenvalues of 
the matrix (|3} as a function of the damping cj> near the crit- 
ical damping. Since the qualitative appearance of graph is 
independent of M for M ^ N, we have avoided putting nu- 
merical values on the axis specific to a particular M. The 
critical damping occurs at cos^ = (1 — sin6)/(l + sin 6) = 
{N - 2^m[n - M))/{N + 2^M{N - M)). 



picts the average number of oracle calls to find the tar- 
get as a function of N and (p in the case M = 1. As 
the damping increases, the Grover rotation shrinks, so 
the number of oracle queries tends to increase. For small 
damping, the average number of oracle queries to find the 
target has the quantum behavior 0{'\fN), while for large 
damping, the number of queries goes like the classical 
expectation 0{N). However, for the critical value of the 
damping, there is a distinct valley that separates these 
two regimes. (The number of oracle calls on the z-axis 
of Fig. 2 is computed assuming the following reasonable 
strategy. Knowing TV and (f>, one plans to execute R calls 
to U. If the external spin flips before the R calls are 
complete, a target has been found with certainty and no 
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FIG. 2: Average number of oracle calls to find the target for 
Af = 1 as a function of A'^ and (p. Note the linear dependence 
on N for large cj>, the \/iV dependence for small 0, and the 
valley separating these two regions at the critical damping. 



FIG. 3: Damping of oscillations of Tr{pZ). They approach 
rather than -1 since p is only the part of the density matrix 
for which the external since has not flipped. 



more calls are needed. Otherwise, the R calls to U are 
followed with one query to the oracle to determine if a 
target has been found. If not, one starts the procedure 
over. By fixing a judicious choice of R, one minimizes 
the expected number of oracle queries. This minimum is 
plotted on the z-axis in Fig. 2.) 

Having pointed out the existance of a critical damping 
value and noted its special character, we now consider 
how the damping can be used in the case of unknown M . 
The qualitative appearance of Fig. 1 is maintained for 
general M <C iV, but the value of the critical damping 
depends upon M through 0. In the absence of knowl- 
edge of Ad , suppose we assume the worst-case scenario of 
fewest targets, M = 1, in which the damping cos0 = 
(1 - sine) /{I + sin 6*) = {N - 2.J M{N ~ M)) / {N + 
2^M{N -M)) = {^/W~^-lf/{^/Tr^+lf is weak- 
est. Fig. 3 shows the progress of the algorithm given 
this choice of 0, for N = 10, 000 with M = 1 and with 
M — 40. The effect of the dissipation is evident; the os- 
cillations into and out of the target state are effectively 
damped even for M substantially greater than 1. The 
value of Tr{pZ) tends to zero, and with high probability 
the external spin flips. Since the external spin signals 
success, it is not necessary to know M to run the search 
effectively. 

Although the choice coscj) = (//V^- l)2/(VlV~^+ 
1)^ is effective for a broad range of M values near M — 1, 
it is not effective when the unknown value of M happens 
to be very large, on the order of N/2. For such large 
values of M, the angle 9 is much greater than this choice 
of 0, and the system approaches the target state long be- 
fore the external spin flips to indicate success. To adapt 
our algorithm so that the spin can flip early when ap- 
propriate, we allow <j) to change from iteration to itera- 
tion. In the first application of U , there is large damping 
= 7r/2 in case M is large. If the external spin has not 
flipped after this iteration of U, this is taken as evidence 
that M must be somewhat smaller, so the damping is 
decreased. For concreteness, for iteration n > 1, we set 



C0S(/)„ = (1 — sin(7r/2n))/(l -|- sin(7r/2n)), which should 
be somewhat near the critical damping for the M satisfy- 
ing R{M) ^ n. Although this choice of 0„ has not been 
optimized, we find that it yields good behavior. 

The results are shown in Fig. 4, which compares the 
average number of iterations before the external spin flips 
to the average number of iterations of the undamped 
Grover's algorithm assuming prior knowledge of M . (The 
average number of iterations of the undamped Grover's 
algorithm is determined as follows. First, one computes 
p(i?), the probability of flnding a target after R calls 
to the oracle given N and M. One then minimizes 
over R the average number of iterations {R + l)/p{R) = 
{R+l)p{R) + 2{R+l)p{R){l-p{R)) + 3{R+l)p{R){l - 
p{R))^ + ... assuming that one performs a verification call 
to the oracle after the R iterations and then repeats the 
whole procedure if the verification turns out negative.) 
In the worst case shown, ignorance of M leads to an ex- 
tra factor of roughly 1.5 oracle calls, although occasion- 
ally damping can actually decrease the number of oracle 
calls. This damping method of coping with ignorance of 
M compares very favorably with other methods such as 
successively applying 2" Grover iterations and measur- 
ing, which behaves quite erractically as a function of M 
sometimes costing an extra factor of 30 or more oracle 
calls, or quantum counting [3], which yields the value of 
M but at a factor far in excess of 1.5 if M is determined 
with reasonable accuracy. 

So far, the damping has been effected "artificially" us- 
ing a single ancilla spin that is subjected to unitary evo- 
lution with the system and then measured. This could be 
the best way to proceed in an actual quantum computer, 
but it is also possible to produce the damping in a more 
traditional fashion by coupling the system to a bath. To 
show this, we first append a flag qubit to the system to 
signal when the target has been reached, taking to 
IV') |0). The flag qubit operators are denoted ax, <Jy, and 
az. Next, we introduce a low-temperature bath, modeled 
as a collection of spins with operators {Sxa, Sy^i, Sz,i}. 
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FIG. 4: Average number of calls to oracle for damped 
search with varying damping COS (pn - — (1 - sin(7r/2n))/(l + 
sin(7r/2n)) as a function of M. We take A'' = 10, 000 and as- 
sume ignorance of M. For each M, the damped search result 
is divided by the average number of calls to the oracle for the 
undamped quantum search assuming that M is known. Note 
that ignorance of M costs at most a factor of 1.5 in oracle 
calls. 



sity matrix for which the flag qubit has not flipped: 
p = - C(l - Z)pI2 - Cp{\ - Z)/2. Here, C 

depends on (j)/9 and on spin correlation functions of the 
bath. This Lindbladt equation leads to dynamics analo- 
gous to those resulting from 

In summary, we have introduced damping into 
Grover's search to mitigate over-rotation past the tar- 
get states when the number of target states is unknown. 
A critical damping value has emerged that divides be- 
tween the classical and quantum regimes. Tuning the 
damping appropriately permits quantum search without 
knowledge of M with only a factor of 1.5 overhead. We 
have presented one promising application of dissipation 
in quantum information, but others can be identified, 
suggesting that this is an exciting avenue for further 
study. 

The author acknowledges helpful discussions with An- 
drew Cross, Daniel Gottesman, Patrick Hayden, and 
Kevin Obenland. 



We can achieve the same results as ^ by revising it to 



U ^ [G 



1 - (T, 1 + cr. 



2 2' 
1 + 5,,, , l-S,,, 



(4) 



The similarity to ([2]) is evident, but there is an ex- 
tra operator [<7x{l + Sz,i)/2 -I- (1 — Sz,i)/2] that flips 
the flag qubit if the bath spin has flipped. We can 
write Q as [/ = exp(-i6ir(l - cr^)/2) J] Jctx^^^ + 
i^]exp(-i(/.5y,,(l - Z){1 - a,)/4) = exp(-i0r(l - 
a,)/2) n^exp(</.(l - Z)/2{a+S+,, ~ a.S.M<^,^±^ + 
"^~^''' ]. Assuming that the bath is cold so that Sz^ — —1 
initially, the last factor can be removed. Now, one can 
think of Grover's algorithm as evolution under a Hamil- 
tonian 0], and since the Grover rotation has the form 
G — exp(— «6'F) it follows that H — Y. To add damping 
to the Hamiltonian, we are motivated by the form of U 
to write H = Y{l-az)/2 + i{(t>/9){l-Z)/2{a+ J2^S+,^ + 
o-~J2iS-,i) + Hbath- Assuming that Hbath allows us 
to make a Markovian approximation, we write a Lind- 
bladt equation [l3| for the portion p of system den- 
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